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We consider a six dimensional brane world model, where the brane is described by a localized so- 
lution to the baby-Skyrme model extending in the extradimensions. The branes have a cosmological 
constant modeled by inflating four dimensional slices and we further consider a bulk cosmological 
constant. We construct solutions numerically and present evidences that the solutions cease to exist 
for large values of the brane cosmological constant in some particular case. Then we study the 
^^ , stability of the model by considering perturbation of the gravitational part (resp. baby skyrmion) 

Cn ■ with fixed matter fields (resp. gravitational background). Our results indicate that the pertur- 

bation equations do not admit localized solutions for certain type of perturbation. The stability 
^ , analysis can be alternatively seen as leading to a particle spectrum; we give mass estimations for 

the baby-Skyrme perturbation and for the graviton. 

in 
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Oh' 
r^,' I. INTRODUCTION 

Theories with extradimensions have been expected to solve the hierarchy problem and cosmological constant prob- 

1^ , lem. Experimentally unobserved extradimensions indicate that the standard model particles and forces are confined 

«ir ' to a 3-brane llMj] . Intensive study has been performed for the Randall-Sundram (RS) brane model in 5 space-time 

r — \ dimensions 0, Q. In this framework, the exponential warp factor in the metric can generate a large hierarchy of 

f~^ ■ scales. This model, however, requires unstable negative tension branes and the fine-tuning between brane tensions 

and bulk cosmological constant. 

There is hope that higher dimensional brane models more than five dimensions could evade those problems appeared 
p^ ' in 5-dimensions. In fact brane theories in 6-dinicnsions show a very distinct feature towards the fine-tuning and 
negative tension brane problems. Warped compactifications are possible in 6 space-time dimensions in the model of 
topological objects such as defects and solitons. In this context abelian strings p-Q were investigated, showing that 
they can realize localization of gravity for negative cosmological constant. (For the positive case, solutions which 
j^ ■ localise gravity are lost [io|-) For the magnetic monopoles, similar compactification was achieved for both positive 
>-j ] and negative cosmological constant |11|. Interestingly, if the brane is modeled in such a field theory language, the 
fine-tuning between bulk and brane parameters required in the case of delta-like branes turns to a tuning of the model 
parameters [9|. 

It is well known that there are two main contexts in which solitons appear in field theories. One is like the strings 
and the magnetic monopoles inabelian and non-abelian gauge theories, and the others are kinds of non-linear type 
models, such as the skyrmions, hopfions [l2|, |l3|. The latter are particularly interesting and have deep insight in their 
non-trivial topological structures. The Skyrme model is known to possess soliton solutions called baby skyrmions in 2- 
dimensional space Il4| -| l 6ll . The warped compactification of the 2-dimensional extra space by such baby skyrmions has 
already been studied [l7[. The authors found regular, static solutions with non-trivial topology which realize warped 
compactification for a negative bulk cosmological constant. Also, a somewhat different model, namely Maxwell gauged 
CP^ type of non- linear a model has been investigated in |18j . 
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Many of previous studies were based on the assumption that the cosmological constant inside the 3-branes is 
tentatively set to be zero; the branes were assumed to be static Minkowski, "flat 3-branes" , despite the fact that our 
universe has a small but positive definite cosmological constant. Thus, addressing the non-zero cosmological constant 
inside the branes has been considered first by Cho and Vilenkin in [19!] and then extended for case of the strings [21| 
and the monopoles [20[. They have studied for both signs of the bulk cosmological constant. In this paper, we 
introduce "inflation" on the baby-skyrmion branes with both signs of bulk cosmological constant. 

Another aim of the present paper is to analyze linear stability of our new solutions by fluctuating all fields. Analysis 
for gravitating thick defects embedded in higher dimensions are found in the literature; for 5-dimcnsions [22|-|2J], and 
for 6-dimensions 0, [25l - [27{ . (Note that the model in [23| is constructed by gravitating multidefects in five dimensions) . 
The studies for thick defects, however, are works in progress since the topological defects used in the literature are 
complicated structures. In this paper, we introduce a general form of perturbation of the metric and the skyrmion 
fields. Employing a simplified ansatz and gauge conditions, we find that two independent equations are sufficient to 
discuss the stablity. We shall solve two of them independently and present typical localized eigenstates of the metric 
and the matter field fluctuations. We shall finally give some speculations for physics of our universe, i.e., SM particle 
spectra, a signature of the CMB and others. 

This paper is organized as follows. In the next section we describe the Einstcin-Skyrmc system in 6-dimensions 
and derive coupled equations for the Skyrmc and gravitational fields. Several types of solutions are found in Sec. III. 
Sec IV is the analysis of the asymptotic behavior of the solutions at the infinity. In Sec V, we present our numerical 
results. We then give a detailed analysis for the gravity and the matter perturbations in Sec. VI. Conclusion and 
discussion are given in Sec. VII. 

II. THE GRAVITATING BABY-SKYRME MODEL IN SIX DIMENSIONS 

The total action for the gravitating baby-Skyrme model is of the form S ~ S'grav + •S'baby- The gravitational part 

5grav=yd'xV^(^^i?-A(6)), (2.1) 

is the generalized Einstein-Hilbcrt gravity action, where A(-g\ is the bulk cosmological constant, and X(e) ~ 87rG(g) = 

On the other hand, the action for the baby-Skyrme model S'baby is given by 

^idMTi) ■ id'^n) ~ ^idMTi X d^nf - KoVin)] . (2.2) 

We will use the convenient notation S'baby = / dPxy/—gCha.hy Here n is a scalar triplet subject to the nonlinear 
constraint n ■ n = 1, and V{n) is the potential term with no derivatives of n. The coefficients K2,4,o in Eq. ()2.2|) are 
the coupling constants in the gravitating baby-Skyrme model. 

A. The Ansatz 

Assuming axial symmetry for the extradimensions, the metric can be written in the following form 

ds^ = M^{p)gl^Jdx''dx'' - dp^ - f{p)d0^ (2.3) 

where p S [0,oo) and 6 £ [0,27r] are the coordinates associated with the extra dimensions. 

We further model a cosmological constant on the brane by considering the following form of the four dimensional 
subspace (described by g^J in Ea. (|2.3p ) 

ds^4) = gl^Jdx^dx" = df ~ S^je^^^'Ux'dx^ (2.4) 

where H{t) is a function of the time coordinate t. 

The ansatz for the scalar triplet n is given by the hedgehog Ansatz [ij] : 

n = {sill f{p) cos{n 9), sin f{p) sui{n 9), cos f{p)). (2.5) 

Let us note that there are some variations j32| for choosing the potential term V{n) in the baby-Skyrme model 
(|2.2|) . Here we use the so-called old baby skyrmions potential, which reads 

y(n) = l-n-n(°°) = l + cos/(p), (2.6) 

where n^°°' = linip^oo ''^■(p, ^) is the vacuum configuration of the baby-Skyrme model. 
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B. Field equations of the model 

In order to rewrite the system in terms of dimensionless quantities, we define 



/^p, L{r):^ !^l{r), n{t) -.^ !^H{t). (2.7) 

K4 V '^4 V '*2 



We also introduce dimensionless parameters ^ according to 
Finally, we introduce 



a:=X(6)'«2, ^ := A(6)K4/k2, fi:=KQKi/Kl. (2.8) 

"(^) ■= 1 + X2(;:y sin' /('^)' ^(^) ■= 1 - JJ(^ sin2 /(r) 

for latter convenience and use the notation /' :~ drf{r). 

After short computations, the baby Skyrme equations and Einstein equations reduce to 

uf" + U— + X + ^ j "/' - (1 + r);^ sin/cos/ + /isin/ = 0, (2.9) 



for the baby Skyrme field and 



3^ + - + 3^ + 3— - 3^ = air, - /3), (2.10) 



A/'2 L" M'L' M" 2n,t.t + 3n\ ^ ^ , ^ 



.M^ M'L' _ n,t,t + 2n^t 

'm2"+ ml ^ AP 



6 — ^+4 3^ ^ ■l=a(Tr-(3), (2.12) 



M" A/'2 n,t,t + 2n\ , ^ , ^ 

'W + 'JP~ ' M^ = "(^« - ^)' ^2-^3^ 

for the Einstein equation. We used the notation H^t '■= dtH{t) and 'H^t,t '■= d^'H{t). 

Let us note that Eq. p.lOp and p. lip are the four-dimensional components of the Einstein equation, while Eq. p.l2p 
and p.lSp are the extra-dimensional components. 

The components of the dimensionless energy-momentum (EM) tensor in Eqs. (l2.10l) - (|2.13l) are given by 

1 1 n^ 

To = -^uf^- -— sinV-M(l+cos/), 

T.-i^/'^-i^sinV-Ml + cos/), (2.14) 

1 1 71^ 

Te - --V /'^ + 9 ;p sin^ / - /i(l + cos /), 

We supplement the system of equations (|2.9p - (j2.13p by the following set of boundary conditions: 

/(0) = -(m-l)^, ./(oo)=^, (2.15) 

where to, n £ Z, for the baby Skyrme field and 

L{0) = 0, L'(0) = 1, M{0) = 1, M'(0) = 0, (2.16) 

for the metric fields. 

The above boundary conditions are required for regularity and finiteness of the energy. 

Note that considering the hedgehog Ansatz (|2.5p under the boundary condition (|2.15|) . one can construct a topo- 
logical charge (or winding number) defined as follows 
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- / n • [dpu X dgn)dpd9 = - [l + (-1)^"'"] £ Z. (2.17) 

t7T J Z 



^ The dimensions of the model parameters are: [A^q)] = M^, [x(6)] = M"' *, [^2] = M'^, [^4] = M'^ [kq] = M^. 



C. Inflating four-dimensional slices 

A solution for 'H(i) is actually trivial. Combining Ea. (|2.1ip and Eq. (|2.10p . leads to 'H^t,t = 0, i.e. H^t — const = "Ho- 
Thus we obtain the simple solution given by 

Hit) ^ Hot (2.18) 

where Ho is a constant; and where we neglected an integration constant. Note that the integration constant can be 
absorbed by a global rescaling of the three dimensional spacial coordinates. 

The result p.lSp follows from the Ansatz for the four-dimensional metric ()2.4p and the important assumption under 
which the gravitational source Chahy is a function of the extra-dimensional coordinates only. 

Therefore we introduce a new dimensionless parameter 

^■.^Hl^^Hl (2.19) 

K2 

and we arrange the components of the Einstein equation in terms of 7: 

3^ + - + 3^ + 3- - 3^ = a(.o - P). (2.20) 

6^-f4 Q-L=a(Tr-fi), (2.21) 

Af2 ML 1/2 V r /^;, y j 



M" M'^ 72 

' h6 6— 

M M2 M2 



4^:r + 6T7J-6T?I = "(^«-/3)- (2-22) 



Note that it is possible to interpret of 7 using the four dimensional effective theory following the lines of ref. [2l| . 
Seen differently, 7 can be interpreted as a positive cosmological constant in the four dimensional subspace of the 
full model, since in this case, g^^ is such that GJ^j, = iH'^g\^^ , where GJ^^ is the Einstein tensor computed with g\^ . 

Note that replacing the four dimensional subspace by another Einstein spacetime satisfying G\jj = iH^g^^j^ , such as 
the Schwarzschild-dc Sitter spacetime, leads to the same equations. 

Another useful quantity is the rescaled Ricci scalar which will be used later and is given by 

L LM M A/2 A/2 ^ ' 

III. SPECIAL SOLUTIONS FOR TOPOLOGICAL VACUUM 

Let us for a moment consider the model with fir) = tt, i.e. where the baby-Skyrme field is a topological vacuum 
configuration. 

It is known that the Einstein equation (|2.20p - (|2.22p for the topological vacuum have several special solutions. Such 
solutions are classified according to the model parameters, mainly /3 and 7. Before reviewing these special vacuum 



solutions, we shall present a useful relation derived in [2l|. From Ea. (j2.2ip . (|2.22p . we have 



M" M'L' ^^ ., . ^ dMjr) 

where Cl is an integration constant. The relation p.ip is quite general except in some special cases, e.g. for A/ = 
and A/' = 0. 

1. Flat branes case: 7 = 

Special solutions for 7 = have been investigated by many authors in Refs. @, [ifl [201 ■ We shortly summarize 
them here. Let Gi, G2 and ro be constants of integrations; for /? = 0, Eqs. ()2.20p - p.22p admit two different types of 
solutions: the string branch 

M^,{t) = C,, Ll{T)^C2{r^ro). (3.2) 



and the Mclvin branch 

M^(r)=Ci(r-rof/5, 



(3.3) 



For /3 > 0, Eqs. ()2.20[) - p.22[) have a set of the periodic solutions given by 



Af°(r)=CiCos^/^/^(r-ro), 



sin,/li^(r-ro) ^3.4) 

L°(r) = C, ^ t-- . 

For /3 < 0, Eqs. (|2.20p - (|2.22l) again admit two different types of solutions; the warped solutions 



M°(r)=Ciexp(±J^r 



L2,M = C2exp(±\/^r 
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(3.5) 



and the divergent solutions 



M°(r) = ClSinh^/^/^(r-^o), 

^ * (3.6) 



-o/.^_ r^„..,u2/5 /-5«/3, 



L:i(r)=C2Cosh^/•^ — ^(r-ro) 



^. Inflating branes case: 7 7^ 

From now on we discuss the inflating branes case, i.e. 77^ 0. We start with a differential equation derived from 
Eas. (|2:20l) - ([2:22l) . given by 

where C is an integration constant. 

Equation p.7p admits two possibilities; Af = and A/' 7^ 0. 

a. The. case of M' = 0. This case is equivalent to M = const = Mc. In this case, equations p.7p and p.ip are 
not valid anymore and we have to consider the full set of Einstein equations. Doing so, we find 



Mir) = ^^, Lir) = ^1 sin U y) ' (3.8) 

which can be understood as follows: the four dimensional slice have a positive curvature, we add extradimensionw with 
positive curvature as well, the curvature of the total spacetime is still positive and related to the positive cosmological 
constant. Note indeed that this solution is well defined only for positive bulk cosmological constant. 

b. The case of M' ^ 0. From Eqs. ()2.20p - p.22p we find a solution of M{r) as a quadrature and a solution of L{r) 
as follows [21I: 



/ A/3 

r — To = / dM 



L(r) = C,— =a 



^ (3.9) 

dM ^ /72A/3-ffA/5 + C 



dr "^ V A/3 



Here we limit the study to the case C = 0, the case with C 7^ involves cllitic functions. For j3 ~ we have a 
cigar- type set of solutions given by 

M,(r)=7(r-ro), L,{r) = Lo = J Cl- (3.10) 

For /? > we again have periodic solutions 
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For /3 < we have diverging solutions 



^ ^ ^ (3.11) 

Lp{r) = LqCOS ( \l -r^ir- tq) 



M.(0 = ^^sinh y^(.-.o) , 

_ ^ (3.12) 

Ld{r) = iocosh ( \/^^(»' - ^0) j ■ 



IV. ASYMPTOTIC SOLUTIONS 

A. Near origin development 

The near origin behaviour of the functions /, L, M subject to the boundary conditions (|2.15p and (|2.16p is given by 

/(r) = -(m-l)7r + /("'(0)r'Vn! + C'(r)"+^ (4.1) 

T(^ , -27~a(/3 + A^-(-l)"/i) r3 ^ ^^ ^ ^, 

Mir) = l + ^l^^^i^±pi^i™^ + 0(r)^ (4.3) 

where /("^ stands for n derivative of /. 

Note that higher order corrections are straightforward to compute. 

B. Large r limit 

In this section, we give the leading asymptotic correction to the functions /, M, L. We focus on the case 7 7^ 0. 
The asymptotic solution is then given by (|3.10p . (|3.12p and (|3.1ip according to the sign of /3. Note however that we 
couldn't find the subleading corrections in the case /3 > 0; in this case, the metric functions are given to the leading 
order in terms of trigonometric functions; it is not possible to neglect such terms. 

1. /3<0 

Considering the topological vacuum solution plus a fluctuation and suppressing subdominant terms in the equations, 
we get 

SM « M-eV-"'9/iO'- + M2"e-V-8"'3/5'-^ (4.5) 

5L « i-eV-^/S/iO'- + ^-g-V-8«/3/5'-, (4.6) 

where 5M and 6f denote the fluctuation around M{r) and f{r) such that f{r) = — (rn — l)7r + Sf{r), M{r) ~ 
Y^— 107/(3a/3) sinh(A/— a/?/10r) + 5M{r)] and where f^ , f2, Mj", M2 ,L^ , L^ are arbitrary constant. 

Clearly, we are looking for the modes with M^ = f2 — L^ = such that 5f and 5A4 are indeed fluctuations. 



Z. /3 = 

Once again, we start from the topological vacuum plus a fluctuation. Here it is possible to solve for the flutcuations 
without further assumptions: 

5/ =: A"x-ii^3(x), (4.7) 

m = Afi° + ^, (4.8) 

8L ^ i? + ^, (4.9) 



where x — T\ Yp- + M; Kn{x) is the modified Bessel function of second kind and /", Af", Afj ,-^1, -^2 ^''"C arbitrary 
constants. In this case, we arc interested to solutions with Li ~ M^ = 0. 

V 5"^ 

Note that for large r, the function / decays as 5 . 

V. NUMERICAL SOLUTIONS 



We solved the system of ordinary differential equations numerically with the solver Colsys |33j for many values of 
the parameters. Due to the large number of parameters, we decided to adopt the following approach: we keep fixed 
the value of the gravitational coupling, the bulk cosmological constant and the winding numbers. Then we vary the 
Hubble factor for different values of the strength of the potential. Note that the model without inflating 4D slices has 



been studied 17[ ; the way we treat the problem allows a direct visualization of the influence of the Hubble parameter 
on the pattern of solutions. 

Before discussing the case 7 ^ 0, we shortly remind the solutions obtained in ref [l4|. Essentially two types of 
solutions were found for any sign of the bulk cosmological constant, for instance, in the case /3 = the branches of 
solutions were discovered, with asymptotic corresponding respectively to flat space (String branch) and one analogue 
to the Melvin universe (Melvin branch). In the following we will present families of solutions for 7 7^ 0, extrapolating 
between these two sets. 

In the analysis, we focused on the surface energy of the skyrmion E and its mean square radius MSR defined as 

/•OO /'OO 

E^2n T^L{r)dr , MSR^l r^f{r)sm^f{r)dr, (5.1) 



where Tg = — tq, see (|2.14p . The mean square radius allows to caracterize the extension of the branc in the transverse 
direction; the more MSR is small, the more the brane is localised. 

Our results are summarized in figure [T] and [2] for a = 0.1, /3 = 0.1 and a ~ 0.05, (3 — 0.1 respectively and 
m ~ n ~ 1. When the value of 7 increases for fixed a, /3, fi the mean square radius decreases while the energy 
decreases for small values of 7 and increases for larger values. In some intermediate values of 7, the energy of the 
skyrmion passes thought a minimum. The minimum occurs at smaller values of 7 when the value of ^ is smaller. 
Note also that increasing values of /i leads to decreasing values of the mean square radius and decreasing values of 
the energy. 

We present typical profiles of the solutions for ?Ti = 2,3 in figure |3] and S] respectively for non vanishing values of 
the parameters. These figures show that there are three possible geometries depending on the sign of the cosmological 
constant (see ref. 12l|): opened (/3 > 0), fiat (/3 — 0) and closed (/3 < 0); all three geometries with angular deficits. 
This effect seems to be a generic feature of a model where the 4 dimensional branes are inflating. 

Pushing the investigation forward, it appears that the inflating baby-skyrmion exists up to a maximal value of the 
parameter 7 (the Hubble parameter) . A second branch of solution exists as well for values of 7 lower than the maximal 
value. The second branch is of the same type as the first one and connects the two types of solutions available in the 
limit 7 — >■ 0. It should be mentioned that we start with the string solution in the case /3 = and increase the value 
of 7 until the second branch is reached, then we decrease 7 along the second branch. Figure [S] and [5] shows some 
relevant numerical parameters caractcrizing the solution (resp. the energy and the square mean radius) with 7 7^ 
and /3 = 0. 

The case n = 2 is however completely different. The solution still exists up to a maximal value of 7, but then the 
solution crashes; in the sense that the second derivative at the origin diverges at the maximal value of 7 (recall that 
the function / behaves like / ^ — (m — l)7r + /2r^ close to the origin). This is shown in figure [7] where the asymptotic 
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FIG. 1: The value of the energy and of the mean square radius as a function of 7 for various values of p and for a = 0.1, /3 = 0.1. 



value of L and the second derivative of / at the origin are shown, we did not find a second branch in this case; if such 
a branch exists, it seems unlikely that it will connect smoothly to the branch we constructed. 

Note that the case n = \ and n = 2 are also different from the geometrical point of view: the scalar curvature 
vanishes at the origin for n = 2 while it goes to a non vanishing constant for n = I. This is illustrated on figure [5] 
and [5] for (3 = 0,n = 1 (resp n = 2) where we show the metric functions, the baby-Skyrme function and the scalar 
curvature; the picture is similar for non vanishing values of 13. 
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FIG. 2: The value of the energy and of the mean square radius as a function of 7 for various values of ^ and for a = 0.05, (3 — 0.1. 
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FIG. 3: The typical profiles of the solutions for different signs of /3 for m — 2. The function / tends quickly to its asymptotic 
value, so it is not possible to distinguish the different profiles for the function /. 
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a=0.0015, 7=0.005,1^=0.1 
m=3, n = l 



FIG. 4: The typical profiles of tfie solutions for different signs of /3 for m — 3. Here again, it is not possible to distinguisli the 
different profiles for the function /. 
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FIG. 5: The value of the derivative at the origin of the baby Skyrme field, the coefHcient Lq of the metric function L in (|3.10p 
for /3 = 0. 
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Gi=0.01, p=0, )j.=0.1 



ci=0.01, p=0, )i=0.1 





m=l, n=l 
FIG. 6: The value of the energy and of the square mean radius for the two branches with /3 = 0, m = 1, n = 1. 




FIG. 7: The asymptotic value of L and the second derivative of / at the origin (in log scale) for jS = Q,m = l,n = 2. The 
second derivative /2 diverges at the maximal value of 7. The window is a zoom if the region close to the maximal value of 7 
showing /2 (not in log scale). 
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FIG. 8: The typical profile of the metric functions, the baby-Skymie field and the reduced scalar curvature for n = l,m — 1. 



y=0.05, a=0.05, n=0.1, p=0 
m=1, n=2 




FIG. 9: The typical profile of the metric functions, the baby-Skyrme field and the reduced scalar curvature for n = 2,m — 1. 
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VI. STABILITY OF THE INFLATING BABY-SKYRME BRANE 



In this section, we consider a general perturbation Hab around the background metric: 

ds^ = A/(r)2 [dt^ - exp{2Ht)dx^) - dr^ - L{rfd9^ + hAB{x^)dx^dx^ , 



(6.1) 



and we choose the transverse traceless gauge V^/i"^^ = h^ = 0. Note that h^^ = g^'^g 
background metric, and the inverse perturbed metric is given to order one in h by g — h 

We wiU focus on the case /3 = 0, already quite involved. 

We parametrize the metric perturbation according to 



AC ^BDj^^^^ g being the 



i-AB 



„nt-iNe 



( M{rfHii{r) H^ir) 

Hi2ir) M{Tfe'"'H22{r) 



H25{r) ': 



H55ir) 



(6.2) 



LirfHeeir) J 



The gauge conditions, \7Ah^^ 

H22 + -^33 + Hii + iJss 



/i^ = reduce to 



1 



Ha 



(6.3) 



L^M^ 



{LM'^L'H^j + L^ {M {4:M'H5j + MH'r,j{r)) - {3H + n)Hij) - iNAI{r)^Hej) = 0, 



1 



jjj^ {LM'^L'iHee - H55) + L^ [M (M'Hn + M' H22 + M'Hs^ + Af'i/44 - 4Af' H55 - MH'^^) 



-{3H + n)Hi5) ^ iNM^H^e) 



0, 



1 



-LML'Hi5 + L^ {M {{3H + n)Hii - H[^ - HH22 - HH33 - HH^i) - ^M' Hi^) 



+iNMHei) = 0, 

where j = 2,3,4,6. From the imaginary part of the gauge conditions, it follows that Hqa = 0, A = 1, . 
then possible to solve the gauge conditions for H55 and i?5„j, m = 1, . . . , 5. 
The matter field is parametrized as follows 

n = (sin(/(r) + 6F{t, r, 6)) sin nO, sin(/(r) + SF{t, r, 9)) cos nO, cos(/(r) + SF{t, r, 9))). 

We parametrize the time and angular dependence according to 

SF{t,r,9)^e^^*-'^^4>{r). 

In this setting the perturbed gravity equation in the transverse traceless gauge are given by 

1 



:(Al) 



ABCD' 



'CD 



- {UgAcgBD + 2Rabcd) h'~^° = a ( 5Tab + - {STqab + THab 



Ah 



AB, 



,6. It is 
(6.4) 
(6.5) 

(6.6) 



where T is the trace of the stress tensor, 6Tab is the first order variation of the stress tensor due to the variation of 
the matter field and of the metric fieldsand where we restored dimensionless functions and parameters. 

We fist work out the gravity equations. The component (2,2), (3,3), (4,4) of the equations (j6.6p are formally the 
same for H22 (rcsp. H33, H44), allowing to set H33 = Hj^y^ = ^22- Furthermore, equations (j6.6p are compatible with 



H12 — -ffi3 — FLu — H: 



23 



H' 



24 



H' 



25 



H: 



34 



H: 



3.5 



0, 



(6.7) 



by looking at the corresponding components of (|6.6p . 

The surviving degrees of freedom are Hii,H22, H^i. It turns out that it is possible to get a single equation out the 
components (1, 1), (2, 2), (5, 5), (5, 1), (6, 6) of (|6.6p by solving these equations for iJ"^, ^22, -^22; -^si- I^i other words, 
we eliminate H^i and it turns out that the equation for H22 and for Hu decouple. Note that solutions to (|6.6p are 
not pure gauge as long as the right hand side does not vanish. 

We will consider the equation for Hu, given by 



h"ir) + Vi{r)h'{r) + V2{r)h{r) + Si{r)(f>{r) + S2{r)(l)'{r) = 0, 



(6.8) 
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V M' 



V2{r) 



n^af'^ sin^ / 1 
2I2 4 



«/ 



/2 






2Af" 
il/ 



(6.9) 
(6.10) 



M2 



3/70; 27 ^2 /J 

-^^ -^ ^ -F2 + ^T " "/^' 



A/2 



M 



2 



A/2 



where we use /i(r) = //ii(r) for shortness, lj = \rT~^ ^^^ where ^i, 5*2 are involved combination of the background 

functions but are straightforward to compute. 

The hnearized equation for the matter fields is quite long and can be obtained from the variation of the reduced 
action, leading to an equation of the following form for the matter field perturbation: 



M^(j)"{r) + Pi{r)(t>'{r) + {P2{r) + m^)(t){r) + Ti{r)h{r) = 0, 



3.11) 



where 



Pi(r):=- 
P2{r) := - 



A/2 



7i2 sin^ / + L2 

A/2 r ,L' 

n2 sin^ / + L2 [ L 



T' 1\T' 

(iV2 „ 3)^2/' sin 2/ + — (n2 sin^ / - L^) - A—in" sin^ / + L^) 
1j i\i 



(6.12) 



/' sin 2/ + n^ cos 2f{{N^ - 3)/'^ + A^^ - l) 



.,M' 



-n2sin2/(/" + 4/'^) - ^iL\N' - I) cos f 



Ti{r) 



where we defined 



g(5V7 + 2c^)(2A/(r)2 + l)/'(r) 
2A/(r)3A/'(r) 

){Zy/^ + uj) (see the following discussion). 



(6.13) 



Note that the system (|6.8p and (j6.1ip constitutes an eigenvalue-like problem, the eigenvalue being essentially given 
by w. 

In fact, m^ appears as the (dimensionless) four-dimensional masses of the scalar four dimensional harmonics ■(/'„: 



V^V^^„ = (92 + i^dt - A)i'ra = -m^A 



(6.14) 



Assuming the particular parametrisation of the perturbations (they don't depend on the four dimensional spatial 
coordinates), m^ indeed reduces to ^uj{3y/j + uj). 
In this case, we easily find the solution of (|6.14p : 



^^^,{-^±mr^>. 



(6.15) 



It follows that the relevant parameter for the stability is m^ ; modes with m^ < leads to tachyonic instabilities while 
modes with m^ > are stable. 



A. Numerical analysis: The formal discussion 



The most obvious boundary conditions are the conditions at infinity, where all the perturbations should vanish: 

0(00) = /i(oo) = 0. (6.16) 

In order not to spoil the topological properties of the baby skyrmion, we need to fix 0(0) = 0. Close to the origin, 
the function h should behave like hK, r^ (the detail of the asymptotic equation is the same as in the next section) . 
In practice, we integrate equations (|6.8p . (|6.1ip in the following way: first, we integrate the equations between and 
some intermediate values, say r^, using a Runge-Kutta algorithm at order 4, with the following boundary conditions: 
0(0) = 0, (/''(O) = (5i, /i'(0) = 0, /i(0) = /iQ for some real values of 5q, Hq, then we integrate the equations backwards 
from a large value of the radial coordinates, say r; to r^, imposing the suitable decay of the functions 0, h and their 
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derivatives, given later. We generate two set of two linearly independent solutions, one between and r„ 

'^^ and h]}, 4)]^^. The general solution is given by 



and hj^ , ( 



y^, the second between r^ and r;, say h]^, yy^ anu /t^, i 



say hi 



L' 



h{r) 



chl- 



dhl, 



for r < rra, 
for r > r™, 

for r < Tm, 
for r > Tm, 



(6.17) 
(6.18) 



where a, b^c,d e M. 
intermediate value r„ 



The solution is smooth if the values of the function and of their derivative match at the 
. The matching condition is expressed by 



det 



hl{rm) hl{r„,) h\.{r,n) hj^{r^) 

'/'i(^m) <l>lirm) (l)R{r,n) 4>i{rm) 

h\'{rm) hl'ir^) h],'ir^) /i| '(r„0 

^l '(rra) 0i '(r-™) ^],^ '{rm) <t'l '{r„,) 



(6.19) 



We computed the decay of the functions /i, 
asymptotic region. In the case /3 = 0, we find 



by solving to the leading order the perturbation equations in the 



^ -[y(W^-l)(^+M)]. 



P^^^ 



Hi-e' 



(6.20) 



where Fi,Hi are normalisation factor. (The detailed analysis will be shown in next subsection.) 

In principle, we are able to integrate equations (j6.8p . (j6.1ip based on the method described above, and actually 
have found some solutions. However, the numerical investigation was plagued by number of difficulties: 

• The coefficients of the system of coupled differential equations are given in terms of numerically computed 
functions. These functions were however computed with a relative precision of order 10~^. 

• Since we integrate the equations from to some maximal value of r, say r^, the boundary condition imposed 
are not exact. We checked that the solution follows the correct decay for larger values of r by integrating the 
solution from re to re + S, S being a real number. It turned out that close to r^ the decay was good, but the 
error due to the cutoff showed up further from re- However, varying re did not influence much the eigenvalues. 

Thus we admit that our results are acceptable as preliminary results which give us an initial guess of the true 
eigenvalues. We will present more details in a further publication and here, we will analyse in detail the stability of 
the matter sector (resp. gravitational sector) in a fixed gravitational (resp. matter field) background in terms of a 
slightly different scheme. 



B. Stability of the Baby-skyrmions with fixed gravitational background 



In this section, we shall solve the eigenequation (|6.1ip where all the gravitational background fields remain unper- 
turbed. This means that terms involving h are eliminated from (j6.1ip . The equation becomes decoupled from h. Since 
the equation involves coefficients to be evaluated numerically, we should of course solve it using numerical methods. 
Before moving to the numerical resolution, we might be able to get some intuitions about sign of the eigenvalues; 
which would give indications on the stability properties of the solutions. 

We follow a scheme employed for the analysis of the six-dimensional Abelian vortex Q. Eq. (|6.1ip can be rewritten 
in the form of a zero-eigenvalue mode Schrodinger equation for the function 



?/)(r) 



exp 






ds 



namely 



where the potential is 



V'" + Vir)^} = , 



V{r) = W{r)' + W^{r) 



Pi 



M2 



M, 



W{r) 



Pi 

2A'P 



(6.21) 



3.22) 



(6.23) 
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Eq. (16.221) can be understood in the context of so called the supersymmetric quantum mechanics 3l| where the function 



W{r) is identified as the superpotential in the SUSY QM. The equation can be rewritten by using the operator 



namely 



A = — + W{r), A^ ^-— + W{r), (6.24) 

ar ar 



^'-t)*^"' ^'^''^ 



Formally, the lowest eigenvalue of the operator AA^ equals to the first excited state of the A^A, obtained by the 
reversing the order of the operator A, A^ . If the ground state of the operator A'^A is a zero mode, i.e.A^Aip ~ 0, one 
easily see that the eigenvalues of the operator AAl' are positive-definite. Note however that the potential in the A^ A 
contains deep negative well, so has many negative eigenvalues. Note also that strictly speaking, Ym is not exactly a 
suspersymmetric potential, due to the presence of the additial (P2 + nn?)/M'^ term. 
Here we plot the effective potential 

VM{r) :== W'{r) + W{rY - ^^ . (6.26) 

In order that (|6.22p exhibits a zero-eigenvalue solution, the potential Vm — m? /M"^ should contains some negatives 
regions. In FigflUl we plot the function Vm(?') for the first few of N with a typical gravitational background. The 
case of m^ > is trivial; by suitably adjusting values of m? we could get the solutions. Those are the eigenvalues 
which might be observed at the LHC or some other probes. On the other hand, for m? < the situation is more 
complicated. If the potential is positive-definite, apparently we have no solution of (|6.22p . thus the possibility of 
tachyonic mode is removed. From this point of view, the solutions for A'' = 0, 1 seem to have no tachyonic mode. 
However, the potentials for N > 2 have a negative pit at the core of the brane, which seems to grow as N increase. 
Thus the tachyonic mode may occurs especially for larger N. For the case of A^ = 2, the potential is shallow but still 
has a negative well; the equation might have a tachyonic mode. In order to go further, however, we have to rely on 
the numerical study. 

The eigenproblem can be solved numerically by the standard predictor- corrector method [30|. The detail of the 
procedure is described in Appendix A. Here we shall examine informations of limiting behavior of the fluctuation (/), 
but first, we will study the near origin and asymptotic behavior of the fluctuation 0. 

1. Asymptotic behaviors 

At the vicinity of the origin, the asymptotic form of solution depends on the winding number n. By using the 
asymptotic solution at the vicinity of the origin (|4.1l) - (|4.3p . one can find the linearized equation of (|6.1ip as 



— + —I 



(6.27) 



which has a solution as the form: 

(j>{r) = c+r^+ + c^r^- , A± = - [(1 - pi) ± ^{1 - pif - Ap^ (6.28) 



where c± are arbitrary integration constants. For n = 1, the coefficients of the equation are 

Pi = --^(2ArV- 57/2-1), p2=A^'-l (6.29) 

where u := /'(O). After a slight examination, one finds the regular solution only the case for N > 2. Similarly, for 
n = 2 

Pi = 1, P2 = -4(7V2 - 1) (6.30) 

so one easily see that for N > 2 the solutions are regular. As a conclusion, the solutions at the origin are regular for 
iV > 2 but not for A^ < 2 . 
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For large values of the radial coordinate, the behavior of the function essentially depends on the sign of /3 := 
KqKh/ K^. For /3 > 0, the metric solutions are periodic. It is an unpleasant feature; the periodic behavior seems to 
interrupt us to find the good asymptotics. Of course it does not necessarily indicate that the solution is unstable but, 
for the time being, we concentrate our analysis on the case /3 < 0. 

Recall that for /3 = 0, ([3A0)) . (|4J)) leads to 



where Lq is a constant. So, sufficiently far from the origin, equation (|6.1ip becomes 



L V L(,/ 7 . 



= 0; 



3.31) 



3.32) 



Since this equation is a kind of Bessel's differential equation, the solutions can be written in terms of Bcssel function 
with real/complex variables. The asymptotic solution is of the form: 







lexp(-J(A.-l)(^+,).) 



(6.33) 



Thus we can find non asymptotically diverging modes of the function cfi for all value of N. (Note that because of the 
condition at the origin, only A^ > 2 is acceptable. ) 

For /3 < 0, from (|3.12p . (|4.4p - (|4.6p . The asymptotic behavior of the solutions is 



7 
/ ~ TT, Af^— sinh^r, L^Lq cosher, ^ := 



-a/3 

To" 



The equation (|6.1ip at the far from the origin is 

sinh^ ^rcj)" + ^(sinh<^r cosher + 4 cosh^ ^r) 






em" 



= 



For sufficient large r, the equation becomes 



(/." + 5C0' - li{N^ 



0, 



the localizing modes are thus 



(r) ^ exp 



h^^Vy +^(^^-1))^ 



(6.34) 

(6.35) 
(6.36) 
(6.37) 



For N > 2 only the minus sign is available while A^ = both signs are acceptable. 



2. The numerical results 



For TV = 2, we have found a tower of positive, discrete eigenvalues for both zero and negative /?. In FigfTl] 
we plot some eigenfunctions cf) of Ea. (|6.11|) for n — 1,2 and /? = or —0.15. The corresponding eigenvalues are 
assembled in Table HI For all cases we have many localized modes; there arc essentially no notable differences for 
the changes of /? or n. On the other hand, for N > 2, the tachyonic modes may exist because of the deep negative 
potential. In Fig ll2l we show the eigenvalues corresponding to the cigcnfunction with no-node for the parameter 
set:a = 0.1, /3 = 0, —0.15, 7 = 0.02, /i = 0.1. At least, within this parameter set, we find that we unavoidably get the 
tachyonic modes except for A^ = 2,3. 

We also study the stability property of the second ( "unstable" ) branch. The eigenvalues are assembled in Table 
im For the second branch solutions, we got the highly localized cigcnfunction for the fluctuation (j>; as a result, the 
eigenvalues are quite higher than the first branches. Note that our obtained eigenvalues are dimensionless one. In 
order to recover the dimensionful one, we need information about the Skyrme parameters K2, K4. In Ref. [17|. we have 
estimated masses of fundamental fermions localized on the warped baby-skyrmion branes with negative cosmological 
constant. In rough speculation, we can extract the value of parameter as \/k2/k4 ~ 10 MeV. From this we find 



^/~n{3H 



Q) = A — m 
K4 



13 TeV : the first branch 



-- 87 TcV : the second branch (6.38) 

for the data of Table [Til Apparently, the first branch is in the TeV scales, but the second branch solution could not 
be observed in experimental facility. 
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C. stability of the metrics with fixed skyrmion bacltground 



If the background skyrniions remain unperturbed, the equation (j6.8p can easily be solved in terms of the method 
which has been described in the case of the baby-skyrmions. This corresponds to omit the terms (/>, 0'. In this case 
we define the eigenvalue l^ := — (cj + ^/j){^ + 2y^). 

1. Asymptotic behaviors 

Similar to the matter field case, at first we examine the asymptotic behaviors of the solutions. By using the 
asymptotic solution at the vicinity of the origin (|4.ip - (|4.3p . one can find the linearized equation of (|6.8p as 



1 N^ 

h" + -h' - -^h = (6.39) 

r r^ 

which has a solution as the form h{r) ^ r^ . Thus, for A^ > 1 we have regular solutions while A^ = arc not adequate 
for the regularity condition at the origin. 

Sufficiently far from the origin, one straightforwardly finds the form of the linearized equation of (j6.8p for /? = 

r^h" + 8rh' + ( 12 + Y72 r^ - ^^Tt'* = 0; (6-40) 



The asymptotic solution is of the form: 



1 / N 



/i(r) - — exp(-— r) (6.41) 



Thus we can find non asymptotically diverging modes of the function h for all value of N . (Note that because of the 
condition at the origin, only A^ > 1 is acceptable. ) 
For /3 < 0, we find the modes 



h{r) ^ exp[(-^c. ± ^^c.2 + ^ ) r] (6.42) 

This solution descreases as r — >■ oo. 

2. The numerical results 

Contrary to the case of the baby-skyrmions, only for N = 2, we have found a tower of positive, discrete eigenvalues 
for both zero and negative /3. For A^ = 1, the solution seems to be tachyonic. 

In FiglTSl we plot some eigenfunctions (f) of Ea. (j6.8p for n = 1,2 and /3 = or —0.15. The corresponding eigenvalues 
are assembled in Table IIIIl Similar to the case of the skyrmions, we plot the effective potential for the gravity in 
FiglH 

1^C(.) := W'{r) + W{rr - ^ . Wir) := ^ (6.43) 

For A^ = 1, the potential has a subtle negatives, so it produces the tachyonic mode. On the other hand, the potential 
of A^ = 2 is positive-definite, so the possibility of tachyonic mode is excluded for the gravity perturbation in this case. 

VII. SUMMARY AND DISCUSSIONS 

The aim of this paper was twofold: first studying the effect of a brane cosmological constant on the six dimensional 
baby-Skyrme brane model and second studying the stability of the corresponding solution. The brane cosmological 
constant was modeled by inflating four dimensional slices, which is relevant in the context of inflationary scenario 
within brancworlds and also at the level of the classical equations: the inflating slices indeed model a four dimensional 
positive cosmological constant. 
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Various types of solutions exist when a bulk cosmological constant is supplemented in the model. However, the 
inclusion of an additional brane cosmological constant has the drastic effect of reducing the number of different types 
of solutions to essentially three, according to the sign of the bulk cosmological constant. The three solutions are 
characterized by specific geometries of the space sustaining the extradimensions (namely opened, flat or closed for a 
negative, null or positive bulk cosmological constant respectively). 

The occurence of these three geometries seems to be indeed a generic feature since the same phenomenon occurs for 
the different model considered in [21|. This idea follows from the fact that these three peculiar solutions are indeed 
vacuum solutions, they are further relevant for the asymptotic form of the solutions in the presence of localized matter 
fields. 

More interestingly, the different types of solutions available in the case of vanishing brane cosmological constant 
seems to be interconnected once the brane cosmological constant is turned on; the two types of solutions available 
in absence of the brane cosmological constant are recovered as limits of these two branches for specific values of the 
parameters. These branches exist up to a maximal value of the brane cosmological constant. 

The second part of this paper deals with the stability of the inflating baby-Skyrme branes constructed in the first 
part. Although the unperturbed equations do not depend on the details of the brane, the stability analysis might be 
different for different four dimensional Einstein spaces. 

We established the equations for the stability of the complete baby-Skyrme branes. Since these coupled equations 
arc particularly involved, we examined in more details the stability issue in two particular simplified situations: the 
one with fixed gravitational backgrounds and perturbed baby-Skyrme field, the second with fixed baby Skyrme field 
and perturbed gravitational field. 

Surprisingly, our results for the case of fixed gravitational background shows that bound states of the baby-Skyrme 
excitation exist only for fluctuations winding more than one time (say N times, recall that there is a cylindrical sym- 
metry in the extradimensions) in the extradimensions. For TV = 2, we found no tachyonic instability while tachyonic 
modes develop for higher values of the winding. The perturbation of the baby Skyrme field can be alternatively seen 
as a particle bound on the brane; we provide a rough estimation of the mass of such particles and it turns out that 
they should be of the order of 10 TeV, i.e. accessible in ongoing experiments. 

The perturbation in the gravitational sector admits bound states for A^ > 1. We found a tachyonic instability for 
N = 1 while the latter disappears for iV = 2. This suggest that the stability issue of the solution depends crucially 
on the winding number of the fluctuation and on the sector explored. We believe that the solutions to the coupled 
problem will confirm the pattern guessed by these preliminary results. 

Let us finally mention some applications of our model: in the context of inflationary models, the particule spectrum 
provided by the stability analysis might have an interesting phenomenological interpretation. As an example, one can 
imagine that the particles emitted by the branes might have a particular signature on the CMB. Also, for unstable 
modes, one could imagine that the instability might turn the inflating brane to non inflating brane, thus providing a 
dynamical end to inflation. Let us note however that these ideas are quite speculative and would deserve deep further 
investigations. 
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Appendix A: Numerical method 

First, we decompose the eigencqauion (|6.11l) into coupled first order equation for u{y),v{y) :— 4>{y) 

dv{y) 



u{y) (Al) 

uiy) + (piiy)^ +P2iy)+P3{y))v{y) (A2) 



dy 

dy V v(y) 
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where 



Pi(.y) 

P2{y) 



y 



y^ 



■n? sin f + L^ 

,,2 



nP sin f + L^ 



ji 

{N"^ - 3)n^f' sin 2/ + —{n^ sin^ f - L^) 
ij 

n^ — f'sm2f + n^ cos2f{{N^ - 3)/'^ + N^ - l) 

1j 



4^in'sin^f + L'] 



-n snil 



.'^I' 



2/(/" + 4/'^)-Mi2(iV2-l)cos/ 



psiy) ■■= -jpv^ 



(A3) 

(A4) 



Here we use a coordinate r :— e^ to become finer the meshpoint at the vicinity of the origin. Essentially the method 
is similar to former described integration method in Sec I VI Al The main difference in this case is the matching 
procedure. We evaluate the matching point ym by following condition: 



Pl{ym.)u{ym.) + {P2{ym) + P3(yrn))w(y,ri) = 



(A5) 



We solve the equation both from the origin and the infinity. In order to match v{y) at y = y™, we multiply a constant 
to the outer solution v{y). Instead of using the Wronski determinant (|6.19p . we introduce an arbitrary (5- functional 
potential at an intermediate value ym'- 



Vs{y) := - 



K(y™)]L'-o ,, , 

d{y-ym) 



v{yni) 



(A6) 



Because, if the (5— functional potential exists, the eigenfunction is continuous at the matching point but its derivative 
is not. Therefore the correction in terms of the first order perturbation 



A^;: 



v*{y)Vs[y)v{y) = ~[v\y,^)t::tlv{y,n) 



(A7) 



efficiently improves the eigenvalue, i.e., the eigenfunction. If the analysis reaches the correct eigenfunction, it no 
longer has discontinuity at all and the computation is successfully terminated. 
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FIG. 10: The potential Vm (fOS)) for n = l with the parameter a = 0.1, P = -0.15, 7 = 0.02, fi = 0.1. 
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FIG. 11: A first few of the fiuctuations (j>- We denote 2S,2P, ••• which indicate N = 2 and the number of nodes. The 
background profiles are of the parameter a = 0.1,7 = 0.02, fi — 0.1. The left is for {m,n) — (1, 1), j3 = 0. The right is for 

(m,n) = (1,2), /3 = -0.15. 
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TABLE I: For the case of Ai" = 2, the eigenvalues m of the equation (|6.1H) for (m, n) 
a = 0.1, 7 = 0.02, fi = 0.1, 13 = 0.0, -0.15. 
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FIG. 12: The lowest eigenvalues m^ (with no-node) for A^ = 2 ~ 10 with the parameter q = 0.1, /3 = 0, —0.15, 7 = 0.02, /i = 0.1. 
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TABLE II: For the case of TV = 2, the eigenvalues m^ of the first and the second branches for (m, n) — (1, l)with the parameter 
a = 0.01, 7 = 0.1, ^ = 0.1, /? = 0.0, -0.1. 
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FIG. 13: A first few of the fluctuations h. We denote 2S, 2P, ■ • ■ which indicate N — 2 and the number of nodes. The 
background profiles are of the parameter a — 0.1,7 — 0.02, /i — 0.1. The left is for (m,n) = (1, 1), /? = 0. The right is for 
{m,n) = (1,2), 13 = -0.15. 
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TABLE III: For the case of TV = 2, the eigenvalues P of the equation (|6.8I) for {m,n) = (1,1), (1,2) with the parameter 
a = 0.1, 7 = 0.02, fi = 0.1, 13 = 0.0, -0.15. 




FIG. 14: The potential Vg estimated by ((03ll for n 
iV = 1,2. 



1 with the parameter a — 0.1, /3 — —0.15,7 = 0.02, /i — 0.1 and for 
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